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Abstract. We compute the spectrum of relic gravitons in a model 
of string cosmology. In the low- and in the high-frequency limits we re- 
produce known results. The full spectrum, however, also displays a series 
of oscillations which could give a characteristic signature at the planned 
LIGO/VIRGO detectors. For special values of the parameters of the 
model the signal reaches its maximum already at frequencies accessible 
to LIGO and VIRGO and it is close to the sensitivity of first generation 

^ experiments. 
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1 Introduction 



In the next few years a number of detectors for gravitational waves, and in particular 
the LIGO and VIRGO interferometers, are expected to start operating in a range 
of frequencies between 10 Hz and 1 kHz. One of the possible signals which could 
be searched, correlating the output of two detectors, is a stochastic background of 
gravitational waves. This background is expected to have different components, with 
different origin: it will get contributions from a large number of unresolved sources 
at modest red-shifts, as well as from radiation of cosmological origin. The latter is 
especially interesting, since it would carry informations about the state of the very 
early universe. 

The basic mechanism of generation of relic gravitational waves in cosmology has 
been discussed in a number of papers, see e.g. refs. [|| ||, the reviews 0, |J and 
references therein. The spectrum can be conveniently expressed using 



Pc dlogf 

where p c is the critical density of the universe, p gw is the energy density in gravitational 
waves and / is the frequency. Particular attention has been paid to the spectrum 
produced in inflationary cosmology. In this case one finds that f2 gw decreases with 
frequency as f~ 2 from 10 -18 Hz to 10~ 16 Hz, and then it is flat up to a maximum 
cutoff frequency / ~ 10 10 Hz. While the frequency dependence of fl gw (f) is fixed, its 
magnitude depends on a parameter of the model, the Hubble constant during inflation. 
An upper bound on the spectrum can be obtained from the measurement of COBE of 
the anisotropy of the microwave background radiation. Via the Sachs- Wolfe effect, a 
large energy density in gravitational waves at wavelengths comparable to the present 
Hubble radius would produce fluctuations in the temperature of the photon cosmic 
background. This gives a limit on || of about 8 ■ 10~ 14 at / ~ 1CT 16 Hz . Since 
for larger frequencies the spectrum predicted by inflation is flat, this bound also holds 
at the frequencies of interest for LIGO and VIRGO. The planned sensitivities of these 
experiments to a stochastic background are of the order of fi gw ~ 5 • 10 -6 , while the 
advanced LIGO project aims at 5 • 10 -11 ||. In any case, the spectrum predicted by 
these inflationary models is too low to be observed. 

Clearly, in order to have a stochastic background which satisfies the COBE bound, 
but still has a chance of being observable at LIGO or VIRGO, the spectrum must grow 
with frequency. A spectrum of this type has been found in ref. in a cosmological 
model suggested by string theory 0, ||, ||. Because of its fast ( ~ / 3 ) growth with 
frequency at low /, the COBE bound is easily evaded, and the most relevant bound 
for this type of spectrum comes from nucleosynthesis. The result is that, for a certain 
range of values of the parameters of the model, the spectrum might be accessible at 
the interferometer experiments, at least at the advanced level, while satisfying the 
existing experimental bounds. 



1 



In ref. || this spectrum has been estimated, using approximate methods, in the 
low- and in the high-frequency limits, and neglecting overall numerical factors. In this 
paper we present a detailed computation of the spectrum, solving exactly the relevant 
differential equations. We fix the numerical factors and we present the frequency 
dependence in the intermediate region. The latter displays an interesting feature: 
it shows a series of oscillations, which might provide a characteristic experimental 
signature. 

As remarked in one must be aware of the fact that it might not be legitimate 
to use field-theoretical methods during the 'stringy phase' of the cosmological model, 
see sect. 2, and large frequencies are indeed sensitive to this phase. However the 
best one can do, at this stage, is to write down a specific cosmological model and see 
what are its predictions. Of course, these predictions should only be considered as 
indicative. 



2 The model 

The low energy string effective action depends on the metric and on the dilaton 
field (we neglect the antisymmetric tensor field). At lowest order in the derivatives 
and in it is given by 

S = ~ 7TT~2 / d4x ^9 \e~* (R + d^dy) - , (1) 



2\l 

where X s is the string lenght and <fi is the dilaton field. The dilaton potential V^(<f>) 
is due to non-perturbative effects and therefore vanishes as exp(— cexp(— 0)) for <ft 
large and negative, with c a positive constant. We consider a homogeneous, isotropic 
and spatially flat background, <fi = 4>(t),ds 2 = dt 2 — a 2 (t)dx. 2 , and we introduce 
conformal time r], dt = a{rj) drj. The pre-big-bang scenario proposed by Gasperini 
and Veneziano f?|, |8|, [| suggests the following choice for the background metric and 
dilaton field. 

For —oo<7]<r] s , with rj s < 0, we have a dilaton-dominated regime with 

at,, = -J- (2) 

H s r] s V a Vs J 

if \ , n V ~ (1 ~ a)Vs /q\ 
nV) = s ~7log • 3 

Vs 

With the values a = 1/(1 + -\/3), 7 = v^3 this is a solution of the equations of motion 
derived from the effective action ([]]) in the absence of external matter [0 . 

At a value r\ — rj s the curvature becomes of the order of the string scale, and the 
lowest order effective action (^ does not give anymore a good description. We are in a 
full 'stringy' regime. One expects that higher order corrections to the effective action 
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tame the growth of the curvature, and both (l/a)da/dt and d<p/dt stay approximately 
constant. In terms of conformal time, this means 



a{rf) 



1 



(4) 



The stringy phase lasts for rj s < r\ < rji < 0. One then expects that at this stage the 
dilaton potential becomes operative and, either with a modification of the classical 
equations of motion due to the dilaton potential [fLOR , or via quantum tunneling |T 



the solution joins the standard radiation dominated solution with constant dilaton, 
which is also a solution of the string equations of motion derived from the action (JJ), 
with external bulk stringy matter |7j. This gives, for rji < rj < rj r , (with i] r > 0) 



a(rj) 



HsVi 



(77 - 2 Vl ) 



(5) 



After that, the standard matter dominated era takes place. We have chosen the 
additive and multiplicative constants in a(r]) in such a way that a{rj) and da/drj (and 
therefore also da/dt) are continuous across the transitions. 

The equation for the Fourier modes of metric tensor perturbations for the two 
physical polarizations in the transverse traceless gauge is fl"2| 



d 2 ip k 
drf 

V(v) 



+ 



a drf 

Inserting the expressions (@-[|) the potential is 

1 



k 2 -V( V )} i; k = 0. 

0/2 ) . 



1 d 2 



v(v) 



(4 ^ _ 1) {7] _ (i _ a ) Vs y 

(V - i)?r 2 , 



-OO < T] < Tj s 
T] s < 7] < 771 
Vl < V < Vr 



(6) 
(7) 



where 2/i = 1 2/3 — 3|, 2v = \2a — 7 + 1|. The exact solutions of eq. (fjD in the three 
regions are 

tpk(v) = \f\v ~ (1 ~ a )Vs\ CHl 2 \k\rj - (1 - a)r) a \) , -00 < 77 < r) s 



Mv) = yf\n\ \A+ H^(k\ V \) + A. HW(k\r)\) 



7] s < 77 < 771 (9) 



B, e tkv - B_ e 



-ikrj 



Vl < V < Vr 
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where H^ 1 ' 2 ' are Hankel's functions. The constants A±, B± can be obtained requiring 
the continuity of the solution and of its derivative. We have chosen the boundary 
conditions so that, at r\ — > — oo, ijj k ~ exp(ikrj). In this case the number of particles 
created per unit cell of the phase space is given by |-B-| 2 . 

Before performing the matching, let us discuss the parameters of the model. The 
two constants a, 7 parametrize the solution in the dilaton dominated phase and there- 
fore they are fixed by the effective action flU): a = 1/(1 + V3), 7 = and then v = 
(anyway, we will write many of our results for generic v). Instead n (or (3) is a free 
parameter which measures the growth of the dilaton during the stringy phase; by 
definition \i > 0. The parameter H s is the Hubble constant during the stringy phase. 
Since in this model the growth of the curvature can only be stopped by the inclusion 
of higher order terms in the string effective action, it is clear that the natural value 
for H s is of the order of the inverse of the string length \ s . If one uses the value 
X 2 S ~ (2/a GUT )Ll l ~ 40 Lpj then the typical value of H s is H s ~ 1/X S ~ 0.15 M pi 
where M v \ is the Planck mass. Finally, there are the two parameters r) s ,r)i. In 
the solution for ipk, and therefore in the spectrum, they appear in the combinations 
fc|77 s |, k\r)i\, where k is the co moving wave number. If we denote by 27r/ the physical 
frequency observed at a detector, we have 2nf = k/a(t pres ), where t pres is the present 
value of cosmic time. Therefore, using eq. @, 

where t eq ~ 3.4 ■ 10 10 4 s is the time of matter-radiation equilibrium, and t pres = 
2/ (3Hq) ~ 2.1- 10 17 tiQ 1 s. The constant h parametrizes the uncertainty in the present 
value of the Hubble constant Hq = 3.2 • 10~ 18 ho Hz, and it cancels in eq. flip]); t± is the 
value of cosmic time when the string phase ends. In this context, the natural choice 
is ti ~ A s . Therefore the parameter rjx can be traded for a parameter fx defined by 

f / U \ / + \ !/ 2 

J J! - A O 1 r,10 T T__ / n S \ I <>1 \ 



*M = ^, . (11) 

The order of magnitude of fi is therefore fixed. Note that at the frequencies of interest 
for LIGO and VIRGO / ranges between 10 Hz and 1kHz, and f / fi is a very small 
quantity. Similarly we can introduce a parameter f s instead of r] s , from k\rj a \ = f/f s - 
This parameter depends on the duration of the string phase, and it is therefore totally 
unknown, even as an order of magnitude. However, since \rji\ < \t) 8 \, we have f 8 <fx- 
To summarize, the model has a dimensionful parameter / s , which can have any 
value in the range < f s < fi and a dimensionless parameter \i > (or equivalently (3 
with 2/i = 1 2/9 — 3|). The dimensionless constants a, v are fixed, a = 1/(1 + V3),v = 
and the dimensionful constants H s , fi are fixed within an uncertainty of about one 
or two orders of magnitude. The constant cb s appearing in eq. (§) drops out from the 
potential, eq. (0), and it is therefore irrelevant for our purposes.^] 

1 For comparison, in ref. g the two parameters which are not fixed are chosen as g s /gi, which in 
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3 The spectrum 



Performing the matching at rj = rj s we get 



.4 



± 



± i 



TcXa\/a 



H^\ax s )H^\x s )-Hf\ax s )H^y 



X s 



J_^_^ H{ 2) { )H{ 1,2) { 



'12) 



where x s = f / f s ; in H^ ,2 \x s ) 1 Hjp(x s ) refers to A + and H^\x s ) refers to A_. 
In deriving these expressions we have used the identity between Hankel functions 
H^\x) Hjp(x) - Hjp'(x) Hl 2 \x) = -4z/(ttx). The constant C appearing in eq. (M) 
has been fixed requiring \ A + \ 2 — \A^\ 2 = 1, which gives \C\ = 1. Next we perform the 
matching at r\ — T]i. However, at the frequencies of interest for LIGO and VIRGO, 
k\r)i\ = f/fi = O(10~ 8 ), and therefore in this second matching we can use the small 
argument limit of the Hankel functions, with a totally neglegible error. This gives 
a relatively simple analytical expression for the coefficient B_ which, apart from an 
irrelevant overall phase, is, for fi ^ 0,0 



7TCt 



2/i- 1 



'J_ 

2f. 



2/i. 



-M-l/2 



Ju 



7 

Js 



(13) 



where J^{z) is the Bessel function. The spectrum of gravitational waves is expressed 
with the quantity 



^gw(/) 



1 dp. 



gw 



-16tt 2 f\B 



Pc dlogf p c 

where p c = 3HqM 2 x / '(87r). Then our result for the spectrum is 



(14) 



where 



^gw(f) = a(n) 



(2tt/ s ) 4 



-hp ( ^ 1 j; U 



fs 



a{n) 



'fx 

Js. 



2/Lt+l 



7 
7. 



5— 2/Li 



2a / 



a 
48 



2 2 ^ (2/i - i) 2 r 2 Gu). 



7 
7- 



(15) 
(16) 



our notations is {f s / fi) 13 , and z s = fi/f s . 

2 For — the small argument limit of the Hankel function is different. The result for /i = 
is the same as eq. ( |l3| ) if one writes 21og///i instead of T(fi) and sets \i — in the remaining 
expression. In the following we write our formulae for /i ^ 0. 
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In the most interesting case v = 
rewrite the spectrum as 

(2vr/ s ) 4 



0, using the identity H { 



(2)'/ 



(2), 



we can 



2^+1 



5-2/i 





+ 



2a f 

Expanding our exact expression for small values of f/f s we get , for v 

(2/.-1) 2 (2vr/ s ) 4 (h 



(17) 



^gw(/) 



192/x 2 « # 2 M 2 V/ s 



(2/xa - 1 + a) 2 + — 

7T Z 



(2/ia — 1 + a) 



'1? 



where 7^ ~ 0.5772 ... is Euler constant. This expression agrees with the result 
obtained in the literature, see eq. (5.7) of ref. 0, apart for the numerical constants 
which cannot be computed using only the approximate solution discussed in refs. || [|. 

Let us observe that the low frequency limit in which eq. flT8| ) holds is actually 
/ -C fs "C fx- If we are interested in the limit / C / s ~ /1 we should not take the 
small argument limit of the Hankel function when performing the matching at 771 . 
Rather, we must keep the exact expression and perform the expansion in the final 
result. If we do not make any assumption on the value of fx/ fs a straigthforward 
computation shows that in the limit f / f s <C 1, f / fx ^ 1 

1 



4/x 



x 



+ M + 



Z7T 



1 — a — 2fia) 1 



111 



1 - a + 2 fiat) 1 



2i / a/ 
2< (, af \ 




-0-1/2 



M-l/2" 



from which we derive again eq. ([PJ) if we now take f s ^fx- As we will see below, the 
graviton spectrum is neglegibly small unless / s -C fx, so the physically relevant limit 
is the one which leads to eq. ([PSD- If we instead consider the spectrum with v > 
a simple calculation gives a low frequency behavior ~ f 3 ~ 2u , without logarithmic 
corrections (the absence of the log / term is due to the different small argument limit 
of Hj) 1 ' 2 ^ for v = and for v > 0.) 

Expanding eq. ( |T7D in the limit / f s (but still / -C fx since eq. ([T7|) holds only 
in this limit) we find instead 



^gw(/) 



4q(/i) (2nf s y f 



2fi+X 



3-2/1 



ir 2 a 



6 



4a(/x) 



7T 2 a HqM^ 




(19) 



which agrees, in the frequency dependence, with the result of ref. ||.[] It is important 
to stress that in the high frequency limit the unknown parameter f s cancels. 

Finally, at sufficiently large /, there is a rather sharp cutoff and the spectrum 
goes to zero exponentially. The cutoff can be obtained computing the spectrum 
without performing the limit / «C /i in the second matching. More simply, the cutoff 
frequency / max can be estimated from fc^ ax ~ |^(?7i)|, which gives 



/max^- VIV-ll/l- (20) 



4 Discussion 



From eq. ( |T9| ) we see that the form of the spectrum depends crucially on whether 
ji < 3/2, /j, = 3/2 or /i > 3/2. Let us consider first the case \i > 3/2. In this case the 
spectrum is a decreasing function of / if / f s . Numerically, 

i^~14.1Cr 6 l ( H ' (21) 

and in eq. (|T9|) this number is multiplied by {fi/f) 2fl ~ 3 ; fi/f is O(10 8 ) at / = 
100 Hz and even larger for smaller frequencies, while 2/x — 3 is positive in this case. 
Therefore, for \x > 3/2, Q gw would violate any experimental bound. More precisely, 
the computation becomes invalid because we should include the back-reaction of the 
produced gravitons on the metric ||. We will therefore consider only < fi < 3/2. 
In this case, the spectrum at low frequencies increases as ~ f 3 log 2 /, and at high 
frequencies is increasing as f 3 ~ 2 ^ (or going to a constant if \i = 3/2.) 

Fig. 1 shows the form of the spectrum for \i = 1.4 and for fi — 3/2. In this figure 
we plot figw(/), measured in units of a(/i) (27r/ s ) 4 /(ifo^pi) {hi fs) 2fl+l , which is 
the overall constant appearing in eq. fllTD, versus f / f s . We see that, compared to 
the low- and high-frequency expansions discussed in 0, the spectrum also displays 
a series of oscillations. Depending on the value of / s , the window available to the 
LIGO and VIRGO interferometers, / between 10 Hz and 1 kHz, may contain many 
oscillations, and this would provide a rather characteristic signature. 

To give an idea of the magnitude of the spectrum, in figs. 2 and 3 we plot h^^f) 
for fi = 1.4 and for \i = 3/2, for a specific value of / s , f s = 100 Hz (this choice of 
parameters is motivated below), and for H s = 0.15M p i and t\ = X s , in the frequency 
range relevant for LIGO and VIRGO. (Note that the quantity of interest for the 

3 For the comparison with ref. B, note that, since 2/i = 1 2/3 — 3|, if 2/3 > 3 the dependence on / 
is - f- 2 ^ = f-W while, if 2/3 < 3, / 3 " 2 ^ = f 213 , which therefore reproduces eq. (3.5) of ref. §. 
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experimentalist is not f2 gw but h^fl^, since h only reflects our uncertainty in the 
quantity which we use to normalize the result.) 

It is also useful to give the result in terms of the quantity h c (f), which is the 
dimensionless strain AL/L produced in the arms of the detector, and is related to 
figw(/) by 



h c (f) ~ 1.3 ■ 1(T 20 V^W) (^7^) • ( 22 ) 

Fig. 4 shows a plot of h c (f) versus / for f s = 10 Hz and for fi = 3/2 and fi = 1.4. 

Let us then discuss what is the best possible result that we can obtain from this 
model, varying the two parameters f s and fi, with < f s < f\ and < fi < 3/2. 
Suppose that we want to detect a signal at a given frequency, say / = 100 Hz. In fig. 5 
we plot f2 gw , from eq. (|ITD , as a function of f s at fixed /. We see that, independently 
of /i, it increases for decreasing f s , and when f s <C / it reaches asymptotically the 
constant value given by eq. ([TP]). So, if we want to detect a signal at a given frequency 
/, the optimal situation is obtained if the value of f s is smaller than /. How much 
smaller is not very important, since as a function of f s , Q gw saturates and practically 
reaches its maximum value as soon as, say, f s < 0.5/. The maximum f2 gw is therefore 
given by eq. (|i~9|), which still depends on the other parameter fi. Since f / fx is a very 
small parameter, we see immediately that the best possible situation is realized when 
/i = 3/2. Note that this means (3 = or (3 = 3, and in the first case not only the 
derivative of the dilaton with respect to cosmic time, but even the dilaton itself stays 
constant during the stringy phase. 

In this case, Q gw reaches a maximum value 



^ max _ 2fc§ (27TAT „ 3 . 10 -r ( H s V ,t ± 
^"gw - 3n F 2 M 2 iU ^0.15 M pl j UJ • 



If fi = 3/2, this maximum value is reached as long as / > f s and therefore, if 
f s is smaller than, say, 10 Hz, it is already reached in the VIRGO/LIGO frequency 
range; after that, the signal oscillates around a constant value (fig. 3). If instead 
fi < 3/2, eq. (|T0f) shows that there is a further suppression factor (///i) 3_2/i , and 
therefore the maximum value is reached only at the cutoff frequency / max ~ fi, 
that is for frequencies around 10 or 100 GHz. However, if fi is not close enough to 
3/2, the suppression factor (///l) 3-2 ^ makes the signal very small at LIGO/VIRGO 
frequencies, unless one uses unnaturally large values of H s ,t 1 (fig. 2). 

Let us then discuss whether this maximum value is compatible with the experi- 
mental constraints mentioned in the Introduction. 

We consider first the nucleosynthesis bound [§, 

J fi gw (/) d(log/) < ^(JV„ - 3) (1 + z^y 1 , (24) 
where 1 + ^eq — 2.32 x 10 4 /ig is the redshift at the time of radiation-matter equilibrium 



and N u is the equivalent number of neutrino species JL5[]. Using the recent analysis 



8 



of ref. [pi , N u < 3.9, we get 

J h 2 n gw (f)d(\ogf)<6.3-10~ 6 . (25) 

For fi = 3/2, figw reaches its maximum value at / ~ / s and stays approximately 
constant untill the cutoff at / ~ fx is reached, and therefore the bound gives 

^- ax log^6.3-l(r 6 , (26) 

Js 

If we take f s ~ 100 Hz (which still satisfies f s < f at the frequency / = 1 kHz 
accessible to LIGO and VIRGO, so that eq. ( p3|) applies), eq. fl26|) gives an upper 
bound on the spectrum predicted by string cosmology 

/^--<3.2-10- 7 , (27) 

which can be obtained from eq. ( p3|) with very natural values of H s ,ti. This number is 
smaller than the planned sensitivity of the LIGO and VIRGO detectors in a first phase 
but well within the planned sensitivity of the advanced project. Smaller values of f s 
give a more stringent bound, but the dependence is only logarithmic: if f s ~ 10~ 7 Hz, 
the maximum value of h^Q gw is 1.5 • 10~ 7 . 
The bound from /xsec pulsars is 

h 2 n gw (f = 10" 8 Hz) < 10~ 8 . (28) 

From fig. 2 we see that, in order to suppress the result at / = 10~ 8 Hz, we must 
have f s ^> 10 _8 Hz, which is well compatible with the condition f s ~ 10 2 Hz required 
before. Even if at / ~ 1 kHz we have h% fi™ w ax = 3.2 ■ 10" 7 , at / < f s we get 



hi Q gw ~ 3.3 ■ 10- J jr log 2 jr . (29) 




We see that a value of, say, f s > 10 Hz is sufficient to bring figw(/ = 10 _ Hz) well 
below the experimental bound. The COBE bound is even more easily satisfied since 
if f s > 10~ 7 Hz the value of f2 gw at / = 10~ 16 Hz is totally neglegible. Fig. 6 shows the 
spectrum for f s = 100 Hz and \x = 3/2 in a large range of frequencies, and compares 
it to the experimental bounds. 

In conclusion, in the most favourable case fi = 3/2 and 10" 7 Hz < f s ~ 100 Hz, the 
relic gravitational waves background predicted by string cosmology at the frequencies 
of LIGO and VIRGO is about hi fi^ x = 3.2 • 10~ 7 , which is smaller than the planned 
sensitivity for coincidence experiments with interferometers in the first phase, but well 
within the sensitivity at which the advanced LIGO project aims. This maximum value 
of f2 gw might also be comparable to the sensitivities which could be reached correlating 



resonant bar detectors such as EXPLORER, NAUTILUS and AURIGA |TJ. At the 
same time, the spectrum satisfies the existing experimental bounds. 
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Figure captions 

Fig. 1 Q gw , measured in units of a(/i)((27r/ s ) 4 /(# 2 M p 2 1 ))(/i// s ) 2 ' t+1 vs. f/f s for fi = 
1.4 and // = 1.5. 

Fig. 2 figw(/) vs. / for // = 1.4, / s = 100 Hz and /i = 4.3 • 10 10 Hz; for comparison we 
also show the low- and high-frequencies limits. 

Fig. 3 Q gw (/) vs. / for /i — 1.5, / s = 100 Hz and /i = 4.3 • 10 10 Hz; for comparison we 
also show the low- and high-frequencies limits. 

Fig. 4 h c (f) vs. / for ii = 1.5 and /i = 1.4, f s = 10 Hz and /i = 4.3 • 10 10 Hz. 

Fig. 5 tt gw vs. f s for /i = 1.4, at fixed / = 100 Hz and f x = 4.3 • 10 10 Hz. 

Fig. 6 figw(/) vs. / for // = 1.5, / s = 10 Hz and /i = 4.3 • 10 10 Hz compared to the 
experimental bounds. 
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